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Abstract. A strong excess in a form of a wide peak in the energy range of 300-800 GeV was 
discovered in the first measurements of the el ectron spectrum in th e energy range from 20 GeV to 
3 TeV by the balloon-borne experiment ATIC (IChang et aUl2008bl) . The experimental data process- 
ing and analysis of the electron spectru m with different crite ria for selection of electrons, completely 



5 independent of the results reported in IChang et al.l (12008b ) is employed in the present paper. The 



Chang et al 



(2008b), but shows that the 



new independent analysis generally confirms the results of < 
spectrum in the region of the excess is represented by a n umber of narrow peaks. The measured 
spectrum is compared to the spectrum of IChang et al.l (I2008bb and to the spectrum of the Fermi/LAT 
experiment . 



1 1 Introduction 



The ATIC (Advanced Thin Ionization Calorimeter) balloon-borne spectrometer was designed to 
measure the energy spectra of elements from H to Fe with an individual resolution of charges in pri- 
mary cosmic rays for energy region from 50 GeV to 100 TeV. The ATIC project had three successful 
flights around the South Pole in 2000-2001 (ATIC-1), 2002-2003 (ATIC-2) and 2007-2008 (ATIC- 
15 4). ATIC-1 was a test flight and is not discussed below. The ATIC is comprised of a fully active 
bismuth germanate (BGO) calorimeter, a carbon target with embedded scintillator hodoscopes, and 
a silicon matrix that is used as a main charge detector. The calorimeter is comprised of 8 layers of 
40 BGO crystals in each layer in ATIC-2 and of 10 layers in A TIC-4. The details of construction 
of the apparatus and the procedures of its calibration are given in 



Guzik et al 



(2004); 



Zatsepin et al 



1 



20 d2004l) : lGanel et alj (l2005l) : lPanov et alj d2008l) . 

It has been shown that ATIC was capable of not only to measure the spectra of cosmic ray nuclear 
components, but also the spectrum of cosmic ray electrons plus positrons (hereinafter referred to as 
the "electron spectrum" for brevity). To separate the electrons from a much higher background of 
protons and other nuclei, the differences in the shower development in the apparatus f or electrons 



25 and nuc lei are used. Spectrum of electrons, measured in this way, was published in IChang et al 



(2008b). The most notable reported detail of the electron spectrum was an "excess" of electrons in 
the range of 300-800 GeV. A possible connection of this "ATIC excess" with the nearby sources of 
cosmic ray electrons such as puls ars and supernova rem nants or annihilation of dark matter particles, 



was discussed for the first time in 
30 discussion in the literature. 
As the result obtained in 



Chang et al . (2008b). Such possibilities provoked a very extensive 



Chang et alj d2008bl) is considered to be a very important, it must be 



tested and confirmed in the ATIC experiment by other methods. The solution of this problem is 
the main objective of this paper. The work has been carried out by the Skobeltsyn Institute of 
Nuclear Physics group of the ATIC collaboration starting from the low-level procedures related to 
35 the apparatus calibration a nd up to the analysis of the final results, completely independently on the 



previous works reported in lChang et al. 



(2008 



ab 



2 Selection of electrons from input cosmic ray particle flux 



The selection of electrons from the input cosmic ray particle flux was performed using the Monte 
Carlo simulations of the showers generated by t he primary electrons in the ATIC apparatus and in the 
residual atmosphere with the FLUKA system dBattistoni et all 120071) . To distinguish the electrons, 
some special quantities are constructed that describe the shape of the shower in the apparatus in 
the longitudinal and transverse directions in such way that these quantities take sharply different 
values for "ty pical" electrons and for "typical" protonsQ. We call these quantities "electron filters". 

Chang et al. ( 2008alb ). where only one filter was used, we constructed five different 



In contrast to 



filters to provide a cross-check of the results and an evaluation of the methodological reliability. The 
following results address to one filter, called Chi, but are confirmed by other filters as well. The basic 
parameters for the construction of the Chi-filter are the relative energy deposits in the layers of the 
calorimeter, Ci —Ei/E (where I = 0,1,2,... is the number of the layer from the top to the bottom, 
and E is the total energy deposit in the calorimeter) and the root mean square widths of the shower 
in the layers i?;. The value of the filter Chi is given by the formula 



Chi = 



\ 



£[(*,- a,)Afl j +£[(*!- f,)K? 



1=0 



1=4 



'The nuclei with z > 2 are easily rejected by the use of the charge detector. 
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where i 7 } = Riy/Ul. The mean values and dispersions of quantities Ri and Fi are calculated for the 
incident electrons by the simulations. The distribution of Chi-values for single-charge particles after 

40 a preliminary selection of events by energy deposits in the layers of calorimeter is shown in Fig. Q] 
The preliminary energy condition consists of a number of cuts for relative energy deposits C; . These 
cuts are the lower limits for Cq , . . . , Cg and the upper limits for C7 , Cg , Cg . They are obtained by the 
simulations in which about 98% of the electrons passed the selection with the applied cuts, while 
about 30% of the protons in ATIC-2 and 70% in ATIC-4 were rejected by these cuts. The energy cuts 

45 are based on the difference of the averaged longitudinal development of the shower in the calorimeter 
for electrons and protons. For the final selection of electrons, we used such cut values for Chi that 
30-25% of the primary electrons were rejected and lost together with the rejected proton background 
(the optimal cut values were different for ATIC-2 and ATIC-4 flights; the efficiency of the selection 
was estimated by the simulations). These losses are taken into account when evaluating an absolute 

50 primary flux of the cosmic ray electrons. 



3 Fine structure in the electron spectrum 

The calorimeter of the ATIC spectrometer is in practice thick for electrons (it is 18 radiation lengths 
for ATIC-2 and 22.5 radiation lengths for ATIC-4), therefore the incident energy of an electron 
can be easily determined by th e total energy dep osit in the calorimeter. The test measurements 



55 on the electron beam at CERN (Ganel et al. 



2005) and the simulations have shown that the ATIC 



spectrometer has a very high energy resolution for the electrons. The resolution is a slow-varying 
function of energy and in the terms of half the line width at half maximum less than 3% at the 
energies of 200-600 GeV. No unfolding procedure is needed to obtain the high-resolution spectrum 
with such a narrow apparatus function since the width of the apparatus function is less than or 
60 comparable with the width of the energy bins for all binnings used. Only an energy-dependent 
scaling factor of about 1.1, evaluated in the simulations, is used to obtain the primary energy of 
the electrons from the energy deposited in the calorimeter. The high value of the energy resolution 
enables us to investigate the electron spectrum for the presence of a structure on the scale of 0.1-0.2 
decade in energy. 

65 It is essential that to detect a structure on the scale of 0. 1 -0.2 decade in energy it is not necessary to 
investigate the "absolute" electron spectrum obtained after the subtraction of the proton background 
(see Fig. [TJ, and accounting for the electron scattering in the residual atmosphere. Neither the 
background nor the scattering of electrons in the atmosphere can lead to a short-period structure in 
the electron spectrum. We have confirmed it by the simulations explicitly: It was shown that the 

70 atmospheric correction produces only an amplitude scaling factor that is slow-varying and almost 
linear function of the logarithm of energy (see Sec. |4] for the further details). On the other hand, 
the protons which can mimic the electrons (Chi ~ 1) and are responsible for the proton background, 



3 



have a very wide energy apparatus function (about 50% of a mean deposited energy), and such an 
apparatus function smears all possible structures (if they are present) in the initial flux of protons. 

75 Consequently, the proton background should form a spectrum without prominent features. 

The spectra of electrons as measured in ATIC-2 and ATIC-4 experiments, without atmospheric 
corrections and the subtraction of proton background, are shown in Fig. [2] in the energy range of 
30-900 GeV with the step of 0.035 decades in energy. It is easy to see a structure in the range of 
200-600 GeV, which is well reproduced in both experiments. Hereinafter, we call this phenomenon 

80 the 'fine structure' . A total statistics per each energy bin of the spectrum is shown in Fig. [3] 

The statistical significance of the observed fine structure is determined by two different factors: 
firstly, by the statistical significance of the presence of non-random structure with the usual ~x 2 - 
criterion for the total ATIC-2 + ATIC-4 spectrum; and secondly, by the statistical significance of the 
correlation (similarity) of structures of the spectra measured separately in ATIC-2 and ATIC-4. The 
evaluation was performed for the region of the spectrum from 200 GeV to 800 GeV. The formulae 
for the calculation of the x 2 -value and a correlation-like function C are respectively: 



where a;,,j/, are integer numbers representing the statistics in the bin number i; Xi,y~i are some 
'smoothed' values of the spectra in the same bin; and erf = y/xl^af = y/yl are the Poisson disper- 
sions for the smoothed spectra; y; denotes the intensity of the sum of ATIC-2 and ATIC-4 spectra 
in the formula for % 2 , and Xi,yi denote the intensities of ATIC-2 and ATIC-4 spectra separately in 

85 the formula for C. Probabilities P x 2 = P(x 2 n d > Xcx P )> Pc = P(C rn d > C cxp ) for the values of x 2 
and C of random spectra to exceed the experimental values Xexp an( ^ Ccxp were calculated. The 
random spectra are the spectra with the Poisson statistics relative to the average values determined 
by smoothed spectra. Quantities 1 — P x 2 , 1 — Pq are statistical significances of the fine structure 
calculated in two different ways. All probabilities are calculated with the Monte Carlo simulations 

90 and are checked approximately against the standard table values where it was possible. 

The problem of calculation of P % 2 , Pc is no way a simple one. One difficulty of the approach is 
that there is some irreducible arbitrariness in the construction of a 'smoothed' spectrum. Therefore, 
one of the main objectives was to investigate the stability of the estimate against an arbitrariness of 
the smoothing procedure. To obtain the smoothed spectra we approximated the primary spectra in 

95 a wide range of energies (from 35 GeV to 1500 GeV) by a cubic spline. To evaluate a stability of 
the statistical significance against the arbitrariness of the smoothing procedure, we studied smooth- 
ing procedures with a very wide range of steps of averaging during the construction of the spline 
(spline steps). We found that there exists a wide plateau in the dependence of an estimated statistical 
significance on the spline step. The estimated statistical significance of the fine structure is almost 
1 00 constant for the spline steps from 0. 12 to 0.24 decades of energy. Such a plateau is possible due to a 
very high amplitude of the fine structure of the spectrum. The details of any 'reasonable' smoothing 
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procedure are not actually important. Our final results are obtained by averaging with spline steps 
0.12, 0.15, 0.18, 0.21, and 0.24. 

The other difficulty is that the estimate is highly dependent on the specific energy binning used in 

1 05 the building of the spectrum. The significance for each spline step (see above) is actually evaluated 
by averaging on a number of different binning with bin widths from 0.015 to 0.035 decades of energy. 
The errors of the final estimates are calculated as standard deviations representing the fluctuations of 
the estimates with different binning and with different 'splining'. The contribution to the error from 
an instability related to the bin size is dominant. 

110 In this way we found that the statistical significance for the correlation is 1 — Pc — (99.69jli 32)%^ 
and for the x 2 -criterion 1 — P X2 = (99.68t°'gg)%. The high statistical significance practically elimi- 
nates a random nature of the observed fine structure. We would like to note that the present estimates 
of the statistical significance should be considered as preliminary ones. It is desirable to develop a 
'binningless' technique to avoid the relatively large statistical errors of the estimates caused by the 

1 1 5 fluctuations of the estimates against binning. This work is currently in progress. 

Several tests were also performed to exclude possible methodological causes of the observed 
structures. The statistics of the proton background in the range of the filter values near the electron 
peak but free of the electron signal (2 < Chi < 3.5, see Fig. [T]) were investigated: no signs of a 
possible structure were found. The different electron filters were studied: all the filters produced 

120 a very similar structure. The spectra for different solid angles and different time periods of the 
experiments were compared: the fine structure was reproduced in all the cases. Thus, no evidence 
that the observed fine structure might be caused by the methodological effect was found. 

If the existence of the fine structure in the electron spectrum is confirmed by independent experi- 
ments, then its most likely source will be the nearby supernova remnants and/or pulsars, but not the 

125 annihilation or decay of the dark matter particles. The idea that the cosmic ray electron spectrum 
in the Te V energy range ma y have features related to a few nearby sources like the pulsars is forty 



years old (Shenetal. 



19701) . The possibility that the deviation from a smooth power-like behavior 



in the spectrum is caus ed by the individual ne arby sources and can be observable in this energy 



range was proposed in 



Nishimura et al 



(1980), in relation to the feasible expe riments. Later on 



130 this is s ue was considered from various viewpoints repea tedly (see, for example. iPohl and Esposito 



(1998); 



Erlvki n and Wolfendald (2002); 



Profumol (120081) ). Specifically, a s t ructur e very similar to 
the one observed in the ATIC experiment is predicted in Ma lyshev et alJ (120091) . where it is es- 
pecially emphasized that such a structure might be used as a signature to distinguish between the 
annihilation or decay of dark matter particles and the other sources of electrons, such as nearby pul- 
1 35 sars. The dark matter c an not be a source of a specific structure represented by several narrow peaks 
(IMalyshev et all l2009l) . The reason is simple and deep: the sources like the pulsars can be treated 
as instantaneous, meanwhile the sources like the dark matter halos or subhalos (clumps) should be 
considered as permanent in time. High energy wing of the spectrum of an instantaneous source can, 



5 



in pri nciple, create a very sha rp peak in the electron spectrum, due to the process of cooling of elec- 



1 40 trons (Malyshev et al 



20091) . therefore a number of sources can produce several peaks as observed 



in the ATIC experiment. On the contrary, the permanent sources like the dark matter clumps would 
mix such peaks with different energies related to the different moments of the time of emission of 



electrons and produce wide distributions in the energy spectrum dKuhlen et al. 



2009) which has little 



in common with the fine structure observed by the ATIC. Even if some 'multi-peak' structure were 
1 45 produced by the dark matter clumps, this structure wou ld be ve ry smooth and sme ared compared to 



the observed ATIC fine structure (see 



Brun et al. 



(2007 FIG.3), 



Cline et al. 



(2010, FIG.9)). 



4 The electron spectrum after proton background subtraction and atmospheric correction 



The spectrum shown in Fig. [2] does not provide a basis for a comparison with the results of other 
experiments, since it does not give a correct absolute intensity of the electron flux. To obtain the 
1 50 correct absolute intensity of the electron spectrum, the proton background must be subtracted, and 
the electron spectrum must be corrected for the scattering in the residual atmosphere (and we neglect 
the secondary atmospheric electron background from a hadronic component of cosmic rays, since it 
is sho wn to be small at the thic kness of a residual atmosphere of about 5 g/cm 2 typical for the ATIC 



155 



flights dNishimura et aUU980l) '). 

The procedure of the subtraction of proton background based on the simulations of proton cas- 
cades in the ATIC apparatus leads to an unsta ble result: the inevi table small errors in the simulations 



lead to large errors in the electron spectrum (Pan ov et al.l j2009). In this paper, we implemented a 



method that is independent of the simulations. The method is based on an approximation of the 
experimental plot of Chi-values for protons (see Fig. Q3 with the simple functions and on the inter- 

1 60 polation of these functions to the coordinates origin. We examined the functions from two different 
three-prpameter assemblages y(x) — Ax 2 exp[— (x / a) a ] and y(x) — Ax/[1 + (x/a) a ), which pro- 
duce a reasonable approximation of the proton peak, but result in somewhat different extrapolations 
of the proton background under the electron peak. The scale of the possible systematic error related 
to the background subtraction is estimated by comparison of the results for different types of func- 

1 65 tions. This method is not rigorous, of course, and should be considered as a qualitative estimate. 
The scale of the proton background is about 5% of the total amplitude of the measured spectrum 
near 40 GeV and about 40% near 700 GeV. The total estimated corridor of possible systematic er- 



rors varies from 



+ 15% 
-16% 



for the energy of 40 GeV to +^ at 700 GeV, and is related mainly to the 



uncertainty in the detection efficiency at low energies, while at high energies it is related to the errors 
1 70 in the subtraction of backgrounds. The systematic errors can not lead to an essential distortion in the 
shape of the spectrum. 

An atmospheric correction is calculated on the basis of the simulations of the scattering of primary 
electrons in the atmosphere using the FLUKA system (Bat tistoni et alj,l2007l) . As suggested in our 
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paper iPanov et al.l (|2009), one does not need to correct the energy at the top of the ATIC apparatus 
1 75 to obtain the energy at the top of the atmosphere since the scattering angles of the secondary gamma 
quanta are very small, and the energy of an electron is recorded in the calorimeter together with 
the energies of almost all secondary gamma quanta. However, these gamma quanta may distort the 
shape of the cascade in the apparatus, which leads to some additional inefficiencies of the filtration 
of the electrons. This was taken into account. The correction factor calculated from the simulations 
1 80 is approximated by a linear function of the logarithm of energy, and varies from 1 .42 at 30 GeV to 
1.26 at 700 GeV, forcing the spectrum to become more steep. Absolute el ectron spectrum mea sured 
in the present paper, along with the ATIC electro n spectrum of the paper IChang et al.l (12008bl) and 



the results of the space spectrometer Fermi/LAT (lAbdo et al.l 120091) are shown in Fig. |4] 



5 Conclusions 

1 85 Our results confirm the existence of the ATIC excess, but this excess is resolved into a fine structure. 
For the energies below 200 GeV, our spectrum is identical in form to the spectrum of Fermi/LAT. 
The fine structure measured by ATIC above 200 GeV can be smoothed by using energy bins as wide 
as in the Fermi/LAT spectrum. Moreover, the ATIC spectrum is related to the part of the South sky 
with declination between —45° and —90°, while the Fermi/LAT spectrum integrates the flux over 

1 90 the whole sky. Taking into account the possible anisotropy of the electron energy spectrum, it would 
be more correct to compare the ATIC spectrum with the spectrum measured by Fermi/LAT for the 
same part of the sky and with narrower energy bins. 
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Chi 



Fig. 1. Black solid line - the distribution of the filter Chi values for single-charge particles (ATIC-2, the energy 
deposit range 100-200 GeV). The narrow peak on the left side is related to electrons, the wide peak is related 
to protons. Blue line represents the cut value to select electron events. The cut value was optimized for the 
best signal to background relation. Red line shows the approximation of the proton distribution for subtracting 
of proton background (see Section [4}. The approximation presented in this picture is from three-parameter 
assemblage y{x) = Ax 2 exp[—(x/a) a ], the region to fit is from Chi = 2 to Chi = 16. 




Fig. 2. The spectrum of electrons without the subtraction of proton background and without an atmospheric 
correction as measured in the ATIC-2 and ATIC-4 experiments. 
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Fig. 3. Total statistics of the electron spectrum in ATIC-2 + ATIC-4 flights in each energy bin. The number of 
events in the range of 200-810 GeV, studied for statistical significance of the structure, is 701. 




Fig. 4. The total ATIC-2 + ATIC-4 electron spectrum of this work aft er subtracting the p roton background 
and after an at mospheric correcti on, the ATIC spectrum of the paper 



Chang et al 



(2008b) and the results 

of Fermi/LA T < Abdo et all 120091) . To simplify the picture, the experimental errors in ATIC spectrum of 



Chan g et al 



1 2008b) are not shown. The estimated corridor of possible systematic errors for the ATIC spectrum 



of the present work varies from 



-Hi*;; 



at 40 GeV to 



+57% 
-46% 



at 700 GeV. 



10 



